Atmospheric escape 


* Importance of escape to planetary evolution 
* Escape mechanisms 
— [hermal escape 


— Non-thermal mechanisms 
* lon pickup 
* Bulk plasma escape along tail 
* Photoelecton-induced ambipolar escape 


e Measurements 


Importance of escapee 


«Mars has no global magnetic field 
«Estimated loss rate ~102 (0.1-0.5 kg) s1 (Lundin et al) — 
significant on solar system timescale 
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Thermal escape mechanism 


* Plot shows the speed 
distributions for oxygen & 
hydrogen atoms in Earth' s 
exosphere, at a typical 
exosphere temperature 
T=1000 K 


* The most probable value is 
vo, Which turns out to be 


2kT 
oars 


* The fraction of oxygen 
atoms above the escape 
velocity (11.2 km s!) is 
negligible, whereas a 
significant fraction of 
hydrogen atoms are above 
the escape velocity. 


* (As the high velocity atoms 

| escape, the distribution 

0 2 4 6 8 10 e 18 quickly returns to this 
Velocity (km s) equilibrium shape, so that 

more atoms can escape.) 


FRACTIONAL POPULATION 


Figure 1.6 Boltzmann distribution of velocities for a molecular ensemble of oxygen atoms 


and hydrogen atoms. Escape velocity v, for earth also indicated. Fig from Salby, Fundamentals of Atmospheric physics 


Thermal (Jeans) escape 


Given a Maxwellian distribution, the upward flux of escaping particles (per unit area) is given 
by the Jeans formula for escape by thermal evaporation: 
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where ) is the number density 


_ JAKT is the most probable velocity, as above 
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-Important factor is the final exponential term. 
* E.g. for Earth, Vese = 11.2 km s-', and assuming exosphere temp T = 1000K: 
— For oxygen atoms, 


Vesc = is the escape velocity 


"" 2x1.38x10 «1000 
i 16x1.66x10 7 


» The exponential term is exp(—11.22/1.022) = exp(—120) = 4 x 103 
which is completely negligible 
— For hydrogen atoms, vy is 4 x larger, giving exp(-7.5) = 5 x 10% . Allows rapid 
escape despite continuous production by H;O photodissociation 


= 1019 m s1= 1.02 km s^! 


Non-thermal escape mechanisms - 1 

* Non thermal mechanisms need to give the escaping particles energies of 
— 0.6 eV/amu on Venus and Earth 
— 0.125 eV/amu on Mars. 


— This is relatively small in comparison with atomic energies and the energies 
which might be gained from an electric field. 


(b) Charge exchange fast neutral escapes 
— Slow neutral + fast ion — fast neutral + slow ion ae 
— By exchanging charge, the VRR 


slow ion 


: i fast ion -- 
fast ion (which was trapped $ trapped 


by the planet's magnetic 
field) becomes a neutral and 
is able to escape. 

— The resulting slow ion is 
trapped by the mag. field. 


(c) Photochemical reactions | 
— These convert the energy of od i 
absorbed X-ray photons to i 
kinetic energy 
* Typical yield is a few eV 


* [his may be enough to 
cause escape 


! neutral 


Non-thermal escape mechanisms - 2 
(d) Sputtering 


— Impact of fast ions knocks atoms out of the atmosphere. The fast particles could 
be solar wind or trapped radiation belt particles 


(e) lon escape 


— lons flow upward along polar field lines and escape down magnetic field lines 


— They gain their energy from ambipolar diffusion, where the interaction of ions & 
electrons in a plasma causes them to diffuse at the same rate 


— The electrons have much higher thermal velocities and can easily escape. This 
sets up an electric field by charge separation and the ions are dragged after the 
electrons 


— Occurs in magnetospheres and also at Mars 
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Non-thermal escape mechanisms - 3 
(f) Electric fields 


— Neutral particles are ionised in the upper ionosphere 

— The resultant ions are pulled out of the ionosphere by electric fields associated 
with auroral activity into the co-rotating plasmasphere above 

— The solar wind sets up a cross-tail electric field (downwards in this diagram) 


— Blobs of cold plasma spin off under the influence of this field as shown and are 
lost into the tail and eventually into the solar wind downstream 


co-rotating blobs of cold plasma 
plasmasphere created by time varying 
solar wind related field 


rotational 
flow 


equatorial slice of magnetosphere 


10.3 Non-thermal escape mechanisms - 4 


f) lon pickup (solar wind scavengin 
— At Mars and Venus the exobase is above the boundary between the solar wind 
flow and the planet. 
— Once the particles become photoionised they can be picked up directly by the 
solar wind flow and carried away. 
— The energy is provided by the electric field in the solar wind. 
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Catastrophic escape mechanisms 


* The following 2 processes are only important in catastrophic 
circumstances involving major impacts: 


(9) Hydrodynamic escape (atmospheric blowoff) 

— When the bulk of a light gas below the exobase moves fast enough it 
streams upward at supersonic speed and drags the rest of the 
atmosphere along 

— This may have occurred due to heating by planetisimal impact during 
the late phase of accretion of the terrestrial planets 


(h) Impacts 
— If the impactor is larger than the atmospheric scale height, the shock 
heated air at the place of the impact will then form a hot bubble that 
rises at more than the escape velocity 


Mars O* lon Outflow 


Lundin et al., 1990 (Phobos data) 
Estimated loss rate ~2x102° s! 


MARS IONOSPHERIC OUTFLOW 
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Fig. 3. Composite diagram of the outflow (n-v,) of 
ionospheric nd 2. 


O* ions from Mars during orbits 1 an 


ASPERA-3 — Preliminary results (IMA) 


Confirmation of the planetary wind - O* and molecular ions 
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Log counts 
per readout 


PEB IMB 


3.0 


X(Rm) 0.21 
Y(Rm) 0.57 
Z(Rm) -0.99 
Alt(Rm) 0.16 


Lundin et al, Science, 2004 
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Figure 3. ZY maps of the integral fluxes (cm? s") of O* , Oj* , CO,*. The direction of the 
interplanetary electric field is show by the red vector. 
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MarsExpress: low rates of loss from high energy ions 
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Q(O+) = 1.610? s1 = 4 g s7 

Q(O,*)21.5:10? s-1 = 8 g s7 

Q(CO,*) = 8:107 s-1 = 6 g s^! , between 30 - 30,000 eV/q. 

Propagated back 3.5 GY, gives 0.2 to 4 mbar of CO, and a few cm of H,O. 


Martian Atmospheric Erosion Rates, S. Barabash, A. Fedorov, R. Lundin, and J-A 
Sauvaud, Science 26 January 2007 315: 501-503 [DOI: 10.1126/science.1134358] 


Escape rates 


Planet Escape rate Reference 

Venus 3dl0^ O*/s * VEX, Fedorov et al., 2008 

Earth 5d10** O*/s (total escape) Geotail, Seki et al., 2001 
43dl0^* O*/s (polar outflow, sol min)  DE-1, Yau etal., 1988 
2240” O*/s (total escape) Cluster, Arvelius et al., 
6740” O*/s (outflow) 2008 
45d 0^ O*/s (inflow) 

Mars (2...9)d10? O*/s ** MEX, Science, 2007 


*  - 10? from Barabash et al., 2008 
** -under revision for low energy escape, e.g. 1025 s! from Lundin et al., 2008 


* [he absolute escape for the magnetized Earth 
may be HIGHER than non-magnetized Mars and 
Venus 


e Strong debate in the community at the moment 
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Fig. 1. Spatial distribution of the proton flux with energy greater than 300eV in cylindrical coordinates: X connects the planet center and the Sun; 
R = vY"+ Z". All scales are in planet radii. The color scale is shown in arbitrary units. Blue color indicates zero flux, and black color shows the absence 
of the measurements. Left panel shows the Martian wake, and right panel shows the Venusian one. Red dashed curve in the left panel shows the theoretical 
position of the magnetosphere boundary after Kallio (1996). Similar dashed line in the nght panel shows the wake boundary at Venus. White dashed 


curves are explained in the text. 
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Fig. 2. The same as Fig. I but for ions of planetary origin with M/Q> 14. 


Fedorov et al., PSS 2009 
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Mars — m/q»14 ion escape controlled by E, Fedorov et al., PSS 2009 
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Mars Express 10 July 2007 — photoelectrons and 
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Fractionof Spectra showing Photoelectron 
Measurements fromELS Sector3 
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Frahm et al, Space Science Reviews 2006 
Frahm et al., Icarus 2010, escape rate of photoelectrons 3.14 x1023 s^ 


The Martian Atmosphere and Water 
Reservoirs 


e Current Atmosphere: 95% CO,, 2.7% N,, 1.6% Ar, 
...., 03 9o H;O. The upper atmosphere also has 
CO, O, NO.... 

e Surface pressure of 6 mb and average 
temperature of 210 K. Below H,O triple point (no 
liquid H,O). 

e Exospheric temperature in 180 - 350 K range; 
exobase is near 200 km. 


e Polar caps have both frozen CO, and H,O. 
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0 180 
Mars Odyssey neutron detection of hydrogen (H,O proxy) at poles 
(Feldman et al., 2004). 


This plus other observations: WEG (Water Equivalent Global Layer) 
of at least 30 m. Atmospheric WEG of 10° m. 


Past Mars - Evolution of Water 


* Evidence for past H,0 amounts far exceeding 
the current amount and in liquid form (unlike 
the present Mars) (c.f. Catling and Leovy, 
2005; Jakosky and Phillips, 2001). 


e Images showing gullies/channels formed by 
flowing liquid (probably about 3.8 Ga). 


* |ssues - what provided sufficient greenhouse 
gas to keep temperature above water triple 
point temp?, and where did all the water go? 
Isotope ratios (particularly D/H ratio of five 
times Earth values) suggests loss of H and O 
(water) over 3 billion years. 


Encyclopedia of the Solar System: Catling & Leovy — Mars Atmosphere: History and Surface 
Interactions 


Fig. 3b. A valley network, centered near 42°S, 92°W. The image is about 200 km across. 
This false color mosaic was constructed from the Viking Mars Digital Image Map (From 
NASA/Lunar and Planetary Institute Contribution No. 1130). 


Current Interpretation (courtesy of B. Jakosky) 


Table 1 Martian isotope ratios and atmospheric loss* 


Isotope ratio Measured valuet Amount lost to space (%)+ 
[From D/H 5 60-74 
Jakosky 98 Ar/ Ar 1.3 ~50-90 
and Poo 1.05-1.07 ~50-90 
Phillips, 'SN/TAN 1.7 ~90 
2002] NS JAG 1.025 ~25-50 


"Values taken from refs 57-59, 62, 77 and 78, and references therein. 

Value estimated, observed or derived for martian atmosphere relative to terrestrial. 
Calculated assuming Rayleigh fractionation. D/H range includes uncertainty in escape 
processes. Other ranges are based on uncertain timing of outgassing relative to escape. 


e Based on current best estimates from SNC meteorites and 
telescopic measurements. 


e Range of amounts lost includes uncertainty in structure of 
upper atmosphere, loss processes, and effects of outgassing 
over time. 


e See also recent review by Lammer et al (Space cience 
Reviews 2008) for Mars and other planets through time 


Mars Evolution contd. 


* Evidence that Mars had a global magnetic field in 
its early history. 

e The dynamo is thought to have turned off about 
3.8 billion yrs ago. Now just remanent crustal B 
remains (MGS magnetometer data). The 
presence of a global field affects atmospheric loss. 

e How did Mars lose its H,O (H and 0) in the past 
and what are the current atmospheric loss 
processes? 


Mars has some remanent crustal magnetism, 
concentrated in the Southern hemisphere -- Mars Global 
Surveyor magnetometer data. 


J.E.P. CONNERNEY ET AL. 


Figure 2. Orthographic projections of the three components of the magnetic field (By, Bg, Bg) at 
a nominal 400 km mapping orbit altitude, viewed from 30 deg S and 180 deg East longitude (after 
Connerney et al., 2001). 


Requirements for Atmospheric Escape 


Atoms/molecules must: 
(1) Be present in the exosphere. 


(2) Have speeds exceeding the escape speed. 


2GM 
ESC R 


Atmospheric Escape Mechanisms 
(1) Hydrodynamic escape (neutral or ion). ** 
(2) Jeans (or thermal) escape (loss of H). 


(3) Nonthermal/photochemical escape (e.g. hot oxygen) 


(4) Nonthermal escape -- planetary ion pickup by the solar wind (or other external 
plasma flow). ** 


(5) Nonthermal escape -- sputtering due to energetic ion impact on an atmosphere. 


LES 


** Solar wind interaction 


Types of Solar Wind Interaction 


e Earth-like (strong intrinsic magnetic field acts as 
obstacle to solar wind flow) 


e Venus-like (ionospheric thermal pressure) 
* Comet-like (ion pick-up and mass-loading) 
e {“real” interactions are a "mixture" of the 


archetypical interaction types - all 3 are relevant 
to Mars} 


Earth-Type Interaction 
The Magnetosphere 
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Figure 1. Schematic of the solar wind interaction with Venus. Interplanetary magnetic : : 
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Pick-up lon Trajectory in the Solar Wind 
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Figure 7.17. Trajectory of an ion in solar wind E and B fields. The ion is initially at rest. 
The ion velocity can be broken down into a gyration part, Ugyro, and an E x B drift part. The 
ion follows a cycloidal trajectory. 


lonization of a neutral by a photon, electron,.. creates a new 
ion which is "picked-up " by electric and magnetic fields and 
partially assimilated into the flow. The mass loaded flow slows. 
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“pickup” by 
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Fig. 1. Illustration of the space environment of a weakly ionized planet [from Luhmann and Kozyra, 1991], showing the 
close-in bow shock, the magnetosheath in which the shocked solar wind plasma flows around the ionopause, and the 
upper atmosphere extending above the ionopause. Ions produced in the solar wind plasma, by processes such as 
photoionization, are "picked up" by the ambient plasma which carries an embedded magnetic field. Because the ions 
gyrate around the field, they may escape into the wake, or they may reenter the atmosphere. 


Ion Pick-up and Loss From Mars Due to the Solar Wind 
From Luhmann and Kozyra (1991); also see Luhmann and Bauer, 
1992) 


Solar Wind Interaction with Mars 


* [he solar wind interaction with Mars has 
characteristics of Venus-like, Earth-like, and 
comet-like interactions although it appears that the 
Venus-like is dominant. 


e Mars has significant crustal magnetic fields in the 
southern hemisphere (but limited global extent). 

* Mars has an extensive exosphere (H and O) 
which leads to pick-up ions and solar wind mass- 
loading. 


Experimental Measurements Related to 
Martian Plasma Environment (key 
spacecraft missions:) 
e Mars-5 
e Viking 1 and 2 
e Phobos-2 
e Mars Global Surveyor 
* Mars Express 
YingHuo-1 (launch 2011) 
* MAVEN (launch in 2013) 


A partially “cometary” interaction partially applies to Mars due 
to its extended hydrogen and oxygen exosphere. 
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Multi-species Global MHD Model: Flow and Magnetic Field 
Results. From Ma et al. (2004). 

Many other global MHD and hybrid code models of solar wind 
interaction with Mars also exist. 
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Figure 4.6. E/q — M/q matrix measured by the ASPERA instrument on March 25, 1989 in the solar 
wind and plasma sheet. The plasma sheet contains mainly ions of the planetary origin. The planetary 
ions gain the energy (E/q) of about 1 keV q^! (from Dubinin et al., 1993a). 


ASPERA experiment on PHOBOS-2 (Dubinin et al., 1993) 
Also see Nagy et al. (2004). Heavy ions in tail (escaped!) 
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Fig. 2. Three examples of energy-m/q spectra for energized ions. Skewed lines indicate nominal 
mass identifications for m/q of 1, 2, 4, and 16, respectively. (A) 27 February 2004, 2- to 7-keV H* 
and ~700-eV heavy ions (O* and O,*) at an altitude of ~290 km. (B) 25 January 2004, ~500-eV 
heavy ions and keV ions with m/q = 2 (possibly He**, H,*, and D+) at an altitude of 900 km. (C) 
22 March 2004, strong energization of heavy ions at an altitude of ~330 km. 


From Lundin et al. (2004)(Science) Several hundred 
eV ions in upper ionosphere - probably escaping. 


Photochemical Escape -- lonospheres of 
Mars (& Venus) 


e hv (ore) + CO, ------ > CO,* + e 
e CO,* + O ------ > O,* + CO 

e OF + CO, ------ > O,* + CO 

e (O5! is the major ion species) 


Dissociative Recombination 
e O;* + e ----- > O + O (source of hot 


* Electron temperature = 1000 - 2000 K and 
ions temperature « 200 - 1000 K 
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Main Non-Thermal Loss Process for Oxygen 


Dissociative Recombination of lonospheric O5* 


O+ + e -> O+ O 


About 50% of the O atoms produced above 
the exobase escape at Mars (Fox and Hac, 1997). 


Prod = 2 a(T,) [O.*] n, 
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Fox and Hac (1997); Fox (2003)- hot oxygen energy distribution 
At several altitudes (200 km exobase). 


Nonthermal/photochemical loss at 
Mars 


(Fox, 1993; Lammer and Bauer, 1991; Fox et al., 1995; Fox and Hac, 1997; 
Kim et al., 1998; Hodges, 2000; Fox, 2003) 


ESCAPE H FLUX = 10° cm? s-' 
(Lyman alpha obsv.) 

(Q = 2 x 1026 5-1) 
ESCAPE O FLUX = 3 x 106 — 2 x 108 cm" s-' 


(Q = 0.3 - 4 x 1026 s^!) 


Summary of Present Oxygen Loss at Mars 


(A) Phobos-2 Taus and Aspera (tail/pick-up): 1-4 
x 1025 51 

(B) Sputtering contribution < 1025s 

(C) Deduced From lonosphere: 2 - 6 x 10495" 

(D) Direct photochemical O escape: 0.3 - 4 x 10% 5-1 

(E) 


E) Global MHD model: = 0.3 x 102°s7 


Discrepancies between models and data (global?) and issue of how 
to extrapolate back in time. “Better” data/models needed — YH-1, MAVEN. 


Table 4 Modeled thermal and non-thermal loss rates of atomic and molecular hydrogen, oxygen, nitrogen and carbon 
species (neutrals and ions) from Mars at present time moderate solar activity conditions (1 EUV), at 2 EUV periods and for 


6 EUV (~ 3.5 Gyr ago) 


Species EUV Jeans Photochem. 
H 1 L5 x 10° 
[1] 
Hh 1 33 x 1024 
[2] 
Ht 1 
Hj 1 
1 28 x 1024 
[4] 
o 2 30x% 
[4] 
z 8.0 x 1075 
[4] 
1 
Ot 2 
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N 1 45 x 1023 
[8] 
oł 1 
C 1 3.0 x 1024 
[10] 
l 80 x 1023 
[11] 
co 2 
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co, 2 
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[5G)] 
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L5 x 1027 
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[3] 
~ 1075 [3] 
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[1] Anderson and Hord (1971), [2] Krasnopolsky and Feldman (2001), [3] Lammer et al. (2003a), [4] Luhmann (1997) for 
I EUV, 2 EUV, 6 EUV also in agreement with Kim et al. (1998) for 1 EUV, [5] Leblanc and Johnson (2002), [6] Penz et al. 
(2004), [7] Lammer et al. (2003b), [8] Fox and Dalgamo (1983), [9] molecular ion outflow is estimated (Carlsson et al 
2006), [10] Nagy et al. (2001), [11] Fox and Bakalian (2001) 
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~1025 4 


Lammer et al., Space Science Rev., 2008 


MAVEN Space Craft images 


Images courtesy of Lockheed Martin 
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MAVEN Will Explore the lonosphere 
and Upper Atmosphere of Mars 
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MAVEN will provide a definitive understanding of the history 
of Martian volatiles. 
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MAVEN (NASA, launch 2013) will measure several aspects of ion escape — but 
low energy escaping population must be measured at higher time resolution 
(~100ms, near escaping ion gyroperiod, to fully study the physics of ion escape) 
— and better electron measurements, multi-point measurements are required 


Summary 

The solar wind interaction with Mars and its atmosphere 
requires the physics represented in all 3 simple 
interaction archetypes - cometary-like (mass-loading), 
Venus-like (ionospheric thermal pressure), and Earth- 
like (intrinsic magnetic field). 

Atmospheric Escape on Mars depends on the details of 
the solar wind interaction as well as on the aeronomy 
of the upper atmosphere and ionosphere. 

YH-1 will make important measurements relevant to this 
problem 

The MAVEN mission will be in an excellent position to 
help unravel escape processes at Mars and to 
provide information for estimating past loss and its 
implications for volatile (e.g., water) inventories on 
this planet. 


Instrument development at MSSL - PNE 
including TechDemoSat (TDS) 7 For 


* Heritage of plasma instrumentation at 
MSSL: Cluster and Double Star 
PEACE, Cassini CAPS Electron 
Spectrometer, Mex and VEx particiption 

e Instrument development for Solar 
Orbiter SWA 


e TDS: UK technology demonstration 
mission 


* Demonstrate miniaturised 
instrumentation under development at 
MSSL 


* Low-resource Electron/lon instrument. 
No mass identification 
* Delivery Oct. 2011, Launch - Q2 2012 


Schematic of TDS instrument 54 
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Enhance auroral activity (T3) 


e Magnetospheric and ionospheric 
disturbances (T1, T3, T4) 
e Disruption of radio communication 
(T1) 
e Disruption of satellite operations 
(T1, T2, T3, T4) 
e Disruption of electrical power grids 
(T3) 
e Malfunction of nuclear power plants 
(T2) 
Timescales (after solar event): 


T1 — electromagnetic radiation reache 


Earth (8 minutes) 


T2 — SEPs reach Earth (20-60 minutes) 
T3 — Plasma reaches Earth (1-3 days) 


space Weather at Earth 
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http://www.swpc.noaa.gov/Media/graphics/Satellite.gif 


T4 — Radiation belt intensifications (T3+2 


days) 


Space Weather at Mars 


e Magnetospheric and ionospheric 
disturbances (T1, T3) 


e Disruption of radio communication 


(T1) iy 

e Disruption of satellite operations + 

(T1, T2, = 

T3) ES 
Timescales: 


T1 - electromagnetic radiation reaches 


Mars (12 minutes) http://mepag.jpl.nasa.gov/science/ 


T2 — SEPs reach Mars (30-80 minutes) 5 Planetary Magnetism/ | 
mars plasmoid Steve Bartlett NASA sm.j 
T3 — Plasma reaches Mars (1.5-4.5 pg 


days) 


